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ABSTRACT 

The clustering pattern in high redshift HI maps is expected to be anisotropic due to two 
distinct reasons, the Alcock-Paczynski effect and the peculiar velocities, both of which 
are sensitive to the cosmological parameters. The signal is also expected to be sensitive 
to the details of the HI distribution at the epoch when the radiation originated. We 
use simple models for the HI distribution at the epoch of reionizaation and the post¬ 
reionization era to investigate exactly what we hope to learn from future observations 
of the anisotropy pattern in HI maps. We find that such observations will probably 
tell us more about the HI distribution than about the background cosmological model. 
Assuming that reionization can be described by spherical, ionized bubbles all of the 
same size with their centers possibly being biased with respect to the dark matter, we 
find that the anisotropy pattern at small angles is expected to have a bump at the 
characteristic angular size of the individual bubbles whereas the large scale anisotropy 
pattern will reflect the size and the bias of the bubbles. The anisotropy also depends on 
the background cosmological parameters, but the dependence is much weaker. Under 
the assumption that the HI in the post-reionization era traces the dark matter with a 
possible bias, we find that changing the bias and changing the background cosmology 
has similar effects on the anisotropy pattern. Combining observations of the anisotropy 
with independent estimates of the bias, possibly from the bi-spectrum, may allow these 
observations to constrain cosmological parameters. 

Key words: cosmology: theory - cosmology: large scale structure of universe - diffuse 
radiation 


1 INTRODUCTION 

Observations of redshifted 21 cm radiation from the large- 
scale distribution of neutral hydrogen (HI) are perceived one 
of the most important future probes of the universe at high 
redshifts, and there currently are several initiatives under¬ 
way towards carrying out such observations. To list a few, 
the GMRT (Swamp et al. 1991) is already functioning at 
several bands in the frequency range 150 to 1420 MHz, while 
LOFAR^, PST^ and SKA® are in different stages of design 
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and/or construction. These observations hold the potential 
of probing the HI evolution through epochs starting from 
the present all the way to 2 ~ 50 when the universe was 
in the Dark Ages (Miralda-Escude 2003). Variations with 
angle and with frequency (or redshift) of the redshifted 21 
cm radiation brightness temperature give three dimensional 
maps of the HI distribution at high 2 (Hogan & Rees 1979). 
The signal is sensitive to the spin temperature Ts of the HI 
hyperfine transition, and at any « the HI will be observed in 
emission or absorption depending on whether Ts is larger or 
smaller than T~, the CMBR temperature. To get a picture of 
the expected signal, we outline below how Ts evolves with 

The evolution of Ts is governed by two tendencies, once 
which drives it towards Try and another towards Tg the HI 
gas kinetic temperature. Ts couples to T-y through a radia- 
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tive transition, and it couples to Tg through atomic collisions 
and through the absorption and re-emission of Ly-a pho¬ 
tons (the Wouthuysen -Field effect, Wouthuysen 1952,Field 
1958 ). The collisional process dominates at high densities 
where it causes Ts to closely follow Tg. The small density 
of free electrons, residual after hydrogen recombination at 
2 ; ~ 1000, transfers energy from the CMBR to the gas and 
maintains Ts — Tg = Ty up to 2 ~ 200. The heat trans¬ 
fer is ineffective at 2 < 200, and the HI cools adiabati- 
cally whereby Ts = Tg < Tj. The collisional process which 
maintains Ts = Tg looses out to the radiative process at 
2 < 30 and Ts approaches Ty. Thus, there is a redshift 
range 200 > 2 > 30 where Ts < Tj and the HI will be seen 
in absorption against the CMBR (Scott & Rees 1990). The 
HI at these epochs largely traces the dark matter, and the 
HI maps would probe the power-spectrum of density per¬ 
turbations at a high level of precision (Loeb & Zaldarriaga 
2004; Bharadwaj & Ali 2004a). The first luminous objects, 
formed at 2 ~ 30, would heat the gas through the emis¬ 
sion of soft X-ray photons (Chen & Miralda-Escude 2004; 
Ricotti & ©striker 2004 )and through Ly-a photons which 
would also couple Ts to Tg through the Wouthuysen -Field 
effect. The HI gas will now be partially heated , and the cou¬ 
pling of Ts to Tg will also be partial depending on the local 
flux of Ly-a photons. This is an additional source of fluctu¬ 
ation in the HI maps, and it may be possible to detect the 
presence of the first luminous objects in the redshift range 
30 < 2 < 15 through this effect (Barkana & Loeb 2004a, 
2004b). The coupling of Ts to Tg is expected to be satu¬ 
rated by 2 ~ 15, and Ts = Tg ^ Tj throughout the gas ie. 
the HI will be seen in emission and the HI maps will again 
trace the dark matter. This changes at 2 ~ 10 when a sub¬ 
stantial fraction of the HI is ionized by the first luminous 
objects. During the epoch of reionization, the HI maps trace 
the size, shape and topology of the ionized regions (Gnedin 
2000; Ciardi et al. 2003; Sokasian et al. 2003a, 2003b; Furlan- 
etto , Zaldarriaga & Hernquist 2004a, 2004b). At low red- 
shifts (2 < 6), the bulk of the HI is in the high column den¬ 
sity clouds seen as DLAs in quasar spectra (Peroux et al. 
2001; Storrie-Lombardi , McMahon & Irwin 1996; Lanzetta 
, Wolfe & syosiTurnshek 1998). These clouds are possibly 
disk galaxies or their progenitors, and they trace the dark 
matter, may be with some bias. The HI will be seen in emis¬ 
sion, and the maps will trace the dark matter (Bharadwaj , 
Nath & Sethi 2001; Bharadwaj & Srikant 2004). 

Though the HI signal in each redshift range will probe 
a different phase of the HI evolution and have it’s own 
distinct signature, there is one thing in common through¬ 
out in that the clustering pattern will be anisotropic. The 
anisotropies in the clustering pattern are the results of two 
distinct effects. The first is the Alcock-Paczynski effect (Al- 
cock & Paczynski 1979, hereafter the AP effect) caused by 
the non Euclidean geometry of space time and the second 
is the effect of peculiar velocities. The AP effect causes ob¬ 
jects which are intrinsically spherical in real space to ap¬ 
pear elongated along the line of sight in redshift space. This 
effect is particularly important at high redshifts (z > .1) 
where it is sensitive to the cosmological parameters which 
it can be used to probe. The two-point correlation function 
of the gravitational clustering of different kinds of objects 
is expected to be statistically isotropic, and this is a natu¬ 
ral choice for applying this test. Galaxy redshift surveys do 


not extend to sufficiently high redshifts for this effect to be 
significant. Quasar redshift surveys extend to much higher 
redshifts and there has been a considerable amount of work 
(eg. Ballinger , Peacock , & Heavens 1996 , Matsubara & 
Suto 1996 , Nakamura , Matsubara , & Suto 1998 , Popowski 
et al. 1998 ,Nair 1999) investigating the possibility of deter¬ 
mining the parameters fim and Ha. The main problem is 
that quasar surveys (eg. SDSS) are very sparse and hence 
they are not optimal for determining the cosmological pa¬ 
rameters (Matsubara & Szalay 2002). Hui & Haiman (2003) 
have used the AP test to constrain cosmological parameters 
using the Lyman-alpha forest. 

The redshift, used to infer radial distances, has a contri¬ 
bution from the line of sight component of the peculiar veloc¬ 
ity and this introduces a preferred direction in the redshift 
space clustering pattern. There are two characteristic effects 
of peculiar velocities. On small scales, the random motions 
in virialized regions causes the clustering pattern to appear 
elongated along the line of sight in redshift space. On large 
scales the coherent infall onto clusters and superclusters, and 
the outflow from voids causes the redshift space clustering 
pattern to appear compressed along the line of sight (the 
Kaiser effect, Kaiser 1987). This effect can be modeled us¬ 
ing linear theory and the anisotropies in the redshift space 
clustering pattern can be used to determine the parameter 
P = /h where Q,m and b are the cosmic mass density pa¬ 
rameter and the linear bias parameter respectively (Kaiser 
1987, Hamilton 1992). The Kaiser effect has been studied in 
the different large galaxy redshift surveys (eg. SDSS , 2df- 
GRS) and a recent investigation of the redshift space dis¬ 
tortions in the 2dFGRS yields a value of /3 = 0.49 ± 0.09 
(Hawkins et al. 2002) at an effective redshift 2 « 0.15 . 

Future redshift survey carried out using redshifted 21 
cm HI radiation will allow the anisotropies arising from both 
the AP effect and the peculiar velocities to be studied to red¬ 
shifts as high as 2 ~ 50 and possibly higher, surpassing the 
potential of any future galaxy or quasar survey. In a recent 
paper Nusser (2004) has studied the AP effect in redshifted 
21 cm maps of the epoch of reionization and he has proposed 
that a determination of the correlation of temperature fluc¬ 
tuations to an accuracy of 20% should allow a successful 
application of the AP test to constrain the background cos¬ 
mological model. 

Barkana & Loeb (2004a) have proposed that measur¬ 
ing the anisotropy of the HI power spectrum arising from 
peculiar velocities provides a means for distinguishing be¬ 
tween the different sources which contribute to HI fluctua¬ 
tions. Further, they propose that it may be possible to sep¬ 
arate the primordial inflationary power spectrum imprinted 
in the dark matter from the various astrophysical sources 
which will also contribute to the HI fluctuations. 

In this paper we re-examine exactly what we hope 
to learn from observations of the anisotropy in HI maps. 
The previous analysis of Barkana & Loeb (2004a) using 
the power spectrum implicitly assumes that the background 
cosmological model is known to a great level of precision 
and hence does not take into account the anisotropies intro¬ 
duced by the geometry (AP effect). We adopt a framework 
which allows high redshift HI observations to be interpreted 
without reference to a background cosmological model. We 
use the HI temperature two-point correlation function which 
deals with directly observeable quantities. Using this frame- 
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work, we quantify the anisotropy and study its dependence 
on the background cosmological model. The AP effect and 
the Kaiser effect differ in their response to variations in 
the background cosmological model. To estimate the rela¬ 
tive contributions from these two effect we also calculate 
the anisotropy ignoring the peculiar velocities. 

Our work focuses on the signal expected from the epoch 
of reionization and the post re-ionization era, these being the 
most promising frequency bands for observations in the near 
future. The HI distribution in these eras is largely unknown, 
and it is expected to differ from the dark matter distribu¬ 
tion. The HI is expected to have a very patchy distribution 
at epoch of reionization. Determining the size, shape and 
distribution of the ionized regions is one of the main forces 
driving the effort towards future HI observations. The post¬ 
reionization HI will probably trace the dark matter with a 
bias. In addition the parameters of the background cosmo¬ 
logical model, the anisotropy in the HI maps is expected to 
also be sensitive to the details of the high redshift HI distri¬ 
bution, the latter being largely unknown. In this paper we 
use simple models for the HI distribution to ask if observa¬ 
tions of the anisotropy tell us more about the background 
cosmological models or the details of the HI fluctuations. 

This paper is organized as follow. In Section 2. we ex¬ 
plain the origin of the anisotropies and how they can be 
quantified. In Section 3. we present the background cosmo¬ 
logical models and models for the HI distribution. Section 4 
contains the results and in Section 5. we present discussion 
and conclusions. 

It may be noted that unless mentioned otherwise, we 
use the values {Q.hh?,h) = (0.02,0.7) throughout. 


2 ORIGIN OF THE ANISOTROPIES 


Observations of the HI radiation at a redshift 2 along the 
direction of the unit vector m will measure 


5Tb{m, z) 


Tb{m, z) - T-y 
TT~z 


( 1 ) 


Tb{z,m.) being the brightness temperature of the HI radi¬ 
ation at the position and epoch when the radiation orig¬ 
inated and Ty the temperature of the CMBR also at the 
same epoch. It is convenient to represent such observations 
as STb{z) where z = zm denotes the space of possible red- 
shifts and directions of observation which we refer to as red¬ 
shift space. Here we assume that the redshift range and re¬ 
gion of the sky under observation are both small. Under 
this assumption, the separation between any two points in 
redshift space is 

Az = Az n -I- z6± (2) 


where n is the line of sight to the center of the region being 
observed, 2 is the mean redshift, A 2 the difference in redshift 
and 6 is the angular separation which is a two dimensional 
vector in the plane of the sky. It may be noted that A 2 n 
and z8i_ are respectively the components of Az parallel and 
perpendicular to the line of sight n, and we introduce the 
notation = Szjz whereby 

Az = 2 (0|| n 4 - 0x) (3) 


We next shift our attention to quantifying the fluctu¬ 
ations in the HI brightness temperature. The temperature 
two point correlation function 


w{9\\ ,9^) = {STb{z) 5Tb{z + Az)) 


(4) 


is the statistical quantity of our choice. Here we have in¬ 
voked the property of statistical isotropy whereby w is in¬ 
dependent of the direction of 9±. The point to note is that 
w(9jj,9±) characterizes the observations solely in terms of 
directly observable quantities, namely angles and redshifts. 
and does not refer to a background cosmological model. It 
is necessary to assume a cosmological model in order to as¬ 
sign a physical position r to a vector z in redshift space. 
Using r( 2 ) to denote the comoving distance corresponding 
to a redshift 2 , 


r(2) 



(5) 


where R( 2 ) is the Hubble parameter, we can assign the vec¬ 
tor r — r(z) m in real space to the vector z = 2 m in red¬ 
shift space. The physical separation Ar (in comoving coor¬ 
dinates) corresponding to the separation Az in eq. (3) is 

Ar = r( 2 ) [p( 2 ) ^ii n -I- 6»x)] (6) 

where p(z) = dlnr( 2 )/dln 2 . The presence of the term p(z) 
in eq. ( 6 ) makes the mapping from redshift space to real 
space anisotropic when p(z) 7 ^ 1. At low redshifts (2 <C 1) 
we have r = cz/H q whereby p{z) = 1 and Ar = (c/Ho) Az 
ie. the mapping from redshift to real space is isotropic. The 
high redshift behaviour of p(z) depends on the cosmological 
model (Figure 1) with the feature p(z) < 1 being common to 
all the models. This introduces an anisotropy, and a sphere 
(Ar)^ = in real space will appear as an ellipsoid (Sy p)^ + 
(9±)^ — (R/r)^ elongated along the line of sight in redshift 
space. This is the Alcock-Paczynski (Alcock & Paczynski 
1979) effect. For a fixed redshift, the elongation depends 
on the cosmological model. For all models, the elongation 
increases with 2 . 

We next discuss how the AP effect will manifest itself in 
ui(0|| ,9±). Assuming for the time being that the fluctuations 
in the HI brightness temperature directly trace the fluctu¬ 
ations in the HI density at the point where the radiation 
originated ie. STb{z) — K 5Hi(r) we have 

w(0||,0x) = A'"eHi(Ar, 2 ) (7) 


where (5hi denotes fluctuations in the HI density, ^Hi(Ar) = 
{&i(r)&i(r + Ar)) is the two-point correlation function of 
5 hi and K is a proportionality factor. We expect the fluctu¬ 
ations in the HI density to be statistically isotropic ie. ^hi 
will not depend on the direction of Ar. It then follows from 
the preceding discussion that w{9^^,9±) will be anisotropic 
when p{z) 7 ^ 1, and the curves of constant w(0||, 9±_) will be 
ellipses elongated along the line of sight. 

The analysis till now has ignored the effect of peculiar 
velocities. The line of sight component of the peculiar ve¬ 
locity makes an additional contribution to the redshift, and 
eq. (5) does not correctly predict the comoving distance. On 
large scales the peculiar velocity field is well described using 
linear perturbation theory which relates it to density fluctu¬ 
ations in the underlying dark matter distribution. The co¬ 
herent flows into over-dense regions and out of under-dense 
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z 

Figure 1. This shows the redshift dependence of the geometrical 
distortion parameter p{z) for different background cosmological 
models. 


ing to the redshift 2 . Here the term phi/ph is the ratio of 
the neutral hydrogen density to the mean hydrogen density. 
Assuming that the hydrogen traces the dark matter and 
that the neutral fraction varies from place to place we have 
PHi{r,z) = ph{z)xhi{z)\1 + S^{v,z)][l + S{v,z)]. Here 5^ 
and 5 refer to fluctuations in the neutral fraction and the 
dark matter respectively, and xhi is the mean neutral frac¬ 
tion. The spin temperature Ts also can vary from place to 
place i .e Ts{yi,z) = Ts{z)[l + 5s(x, z)]. The last term in 
eq. (10) incorporates the effects of linear peculiar velocities 
where v{r, z) is the line of sight component of the peculiar 
velocity and dv/dr is its radial derivative. We then have, at 
linear order, the fluctuating part of ?7Hi(r, z) to be 


?fHi(r,2) 


Xhi 



[5(r,2) -|-5s:(r,z) 


(1 -h z) dv{r, z) 
H{z) dr 


+ p^5s(r,2) 


( 11 ) 



Figure 2. This shows the redshift dependence of the peculiar 
velocity distortion parameter /(z) for different background cos¬ 
mological models. It may be noted that for the SCDM model 
/(z) = 1 throughout. 


regions leads to the compression of structures along the line 
of sight (Kaiser 1987). The introduces an an anisotropy pat¬ 
tern which is opposite in nature to the AP effect. 

Incorporating the effect of the coherent flows to lin¬ 
ear order, the fluctuations is the brightness temperature 
of the HI radiation has been calculated (Bharadwaj & Ali 
2004a,Bharadwaj & Ali 2004b) to be 


is the mean brightness temperature of redshifted 21 cm ra¬ 
diation which depends only on z and the background cos¬ 
mological parameters , and 


77Hi(r 


Phi 

PH 


(‘-1) 


1 - 


{I + z) dv 


H{z) dr 


( 10 ) 


is the “21 cm radiation efficiency in redshift space”. The 
term pHi(r, z) should be evaluated at r, the position in real 
space corresponding to z - calculated using eq. (5)) without 
the effect of peculiar velocities, at the epoch correspond- 


This is most conveniently expressed in Fourier space as 

Pk 


77hi(x, z) = 
where 


(27r)' 


■pHi(k, 2 )e 


( 12 ) 


T-y 


PHi(k,z) = a:Hi( 2 )<j (1 - ^) [A(k, 2 )-I-A,„(k, 2 ) 


+ /P^A(k,2)] -b ^A^kz) 

J- S 


(13) 


Here p = n ■ k/A: is the cosine of the angle between the 
Fourier mode k and the line of sight n towards the center of 
the part of the sky under observation, and / = dlnD/dlna 
where D{a) is the growing mode of density perturbations 
(Peebles 1993) and a is the scale factor. The term involving 
/ arises from the peculiar velocities. 

The HI correlation function ^Hi(Ar, z) = 
(pHi(r, z)» 7 Hi(r-b Ar, z)) and the HI power spec¬ 
trum PHi(k, z) defined as (? 7 Hi(k, 2 )pHi(k , z)) = 

(27r)®5®(k — k')PHi(k, z) are related as 

eHi{Ar,z)= I (14) 

where the HI power spectrum is 

PHr(k,2) = xliiz)l^(l-^y [(l + fp^fPs-s{k,z) 


5Ti,{z) = T(z) X pHi(r, z) 


(8) 

-b 

where 



-b 

T(.-) + 

/0.7\ Ho 

V h H(z) 

(9) 

-b 


B[(l + /p")P«-s(,«) 

J- S 


P.- 


4-1) 

(k,z)]+ (E)' P,_4fc,z) 


(15) 


Here Ps-s{k, z), P^-j;{k,z) and Ps-s{k,z) are the power 
spectra of the three different sources which contribute to 
the HI fluctuations, the dark matter, neutral fraction and 
spin temperature fluctuations respectively. The cross power 
spectra Ps-j;{k, z), Ps-s{k, z) and Px-s{k, z) quantify cross¬ 
correlations, if any, between these different sources. The 
presence of peculiar velocities makes the HI power spectrum 
PHi(k, z) anisotropic through the terms involving /p^ which 
appear in eq. (15). It should be noted that /p^ appears only 
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with the sources of HI fluctuation that are correlated with 
the dark matter fluctuations, and the contribution to the 
HI power spectrum from components not correlated with 
the dark matter are isotropic. 

The anisotropy of PHi(k, 2 ) is reflected in ^Hi(Ar, 2 ) 
which depends on how Ar is aligned with respect to n. 
Further, ^Hi(Ar, 2 ) depends only on the even powers of 
n • Ar/Ar. Shifting our focus back to w(6||, 6±) we have 

«;(0||,0x) = r'W.CHi(Ar,2) (16) 

The correlation in the HI brightness temperature fluc¬ 
tuations, w(0||,0x) is anisotropic due to two reasons, first 
because the mapping from to Ar is anisotropic (the 

AP effect) and second because ^hi is intrinsically anisotropic 
due to the effect of peculiar velocities. Here we investigate 
whether it will be possible to observationally distinguish 
these two distinct sources of anisotropy. Both these effects 
depend on the background cosmological model, and we in¬ 
vestigate to what extent observations of the anisotropy can 
be used to determine the cosmological parameters. Finally, 
the anisotropy in the HI power spectrum depends on the 
properties of the different sources which contribute to HI 
fluctuations. We investigate to what extent observations of 
anisotropies in the HI can disentangle the contributions from 
the different sources. 


3 MODELS 

3.1 Background Cosmology 

We have considered a few representative cosmological mod¬ 
els, all of which are spatially flat and have two components 
namely dark matter and dark energy. The possibility that 
the dark energy equation of state evolves with time has, of 
late, received a considerable amount of attention (eg. Pee¬ 
bles & Ratra 1988 , Riess et al. 1998 , Perlmutter et al. 
1999 , Sahni & Starobinsky 2000 , Sahni & Wang 2000 
, Saini et al. 2000, Carroll 2001, Huterer & Turner 2001, 
Perlmutter 2003, Linder 2003, Alam et al. 2004 ) . Here we 
consider a possible parameterization of the variable dark en¬ 
ergy (VDE) model introduced by Wang & Tegmark (2004). 
The Hubble parameter is of the form 

H{z) = Ho [H™o (1 + zf -1(1- H™o)/db( 2 )]^/". (17) 

where Hq, Hmo, and 1 — QrnO are the present values of 
the Hubble parameter, the dark matter density parame¬ 
ter and the dark energy density parameter respectively, and 
IdeIz) = Pdb(«)/Pdb(0) = /oo + (1 - /oo) (1 + 2 )®^^+“*^ is 

the dimensionless dark energy density function. This VDE 
model has two extra parameters foo and Wi compared to 
the usual LCDM model, and the parameters are restricted 
to the range Wi > —2 and foo > 0. The VDE model re¬ 
duces to to LCDM model if either foo = 1 or Wi = — 1. In 
addition to the SCDM model (HmO = 1) and the LCDM 
model with flmo = 0.3, we have considered the VDE model 
with Q,m 0 = 0.4 and {foo,Wi) = (1.2,—1.8). In the VDE 
models which we have considered, foEiz) increases from 
foEiO) = 1 and saturates at foo = 1.2 around 2 ~ 3. The 
value of foE reaches the saturation value at a lower red- 
shift is Wi is decreased. Figures 1 and 2 show the behavior 
of p(z) and f{z) respectively for the models considered here. 


The points to note are (1.) Though in all models p{z) falls 
rapidly up to 2 ~ 5 beyond which it decreases gradually, 
the values of p{z) are different in each of the models, and 
(2.) The values of f{z) differ from model to model only for 
2 < 4, and f{z) ~ 1 at larger redshifts where the universe is 
dark matter dominated. 


3.2 The HI distribution 

We focus our attention on two different epochs which we 
discuss below. 


3.2.1 Reionization 


We adopt a simple model for the hydrogen distribution 
where there are spherical ionized bubbles of comoving ra¬ 
dius R and the region outside the bubbles is completely 
neutral. The total hydrogen density is assumed to trace 
the dark matter distribution. The reionization is believed 
to have been caused by the first luminous objects which are 
expected to be highly clustered. We incorporate this by as¬ 
suming that the centers of the ionized bubbles are biased 
with respect to the dark matter distribution with a bias be. 
Further, we assume that the hydrogen has been heated be¬ 
fore reionization ie. Ts 3> Ty and consequently the HI will 
be observed in emission. The fluctuations in the neutral frac¬ 
tion (5x) has two parts, one which is correlated with the dark 
matter fluctuations and another, arising from the discrete 
nature of the ionized regions, which is uncorrelated with the 
dark matter distribution. The spin temperature fluctuations 
{Ss) do not contribute when Ts Tj. The power spectrum 
for the resulting HI distribution is (Bharadwaj & Ali 2004b) 


PHi(k, 2) = ^xm(l + frd) 

^ flW\kR) 
nm 


befvW(kR)]' Ps_s(k,z) 

(18) 


where W{y) = (3/y®)[sin(i/) — ycos{y)] is the spherical top 
hat window function, fv is the fraction of volume which is 
ionized, hui is the comoving number density of the spheres 
and /v = (1 — *Hi) = ^TvR^nni. 

of spherical HII region The first term which contains 
Ps-s{k) arises from the clustering of the hydrogen and the 
clustering of the centers of the ionized spheres. The second 
term which has 1 /hui is the Poisson noise due to the discrete 
nature of the ionized regions. The latter is not correlated 
with the dark matter. 

This model has a limitation that the HI density is neg¬ 
ative in a fraction ^ fy /2 of the total volume where ionized 
spheres overlap. The possibility of the spheres overlapping 
increases if they are highly clustered. This restricts the range 
of fv and bo where this model is meaningful. 


3.2.2 Post-reionization 

At low redshifts the bulk of the HI is in the high column 
density clouds which produce the damped Lyman-a absorp¬ 
tion lines observed in quasar spectra (Peroux et al. 2001 
, Storrie-Lombardi , McMahon & Irwin 1996 , Lanzetta , 
Wolfe & syosiTurnshek 1998). These observations currently 
indicate Qgas{z) , the comoving density of neutral gas ex¬ 
pressed as a fraction of the present critical density, to be 
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nearly constant at a value Q,gaa (z) ^ 10“® for z > 1 (Peroux 
et al. 2001). The damped Lyman-aZp/ia clouds are believed 
to be associated with galaxies which represent highly non¬ 
linear overdensities. It is now generally accepted from the 
study of the large scale structures in redshift surveys and N- 
body simulations that the galaxies (or nonlinear structures) 
are a biased tracer of the underlying dark matter distribu¬ 
tion (Kaiser 1984 ;Dekel & Lahav 1999;Mo & White 1996 
;Taruya & Suto 2001 and Yoshikawa et al. 2001). On the 
large scales of interest here it is reasonable to assume that 
these HI clouds trace the dark matter with a constant, linear 
bias b. 

Converting i}gas to the mean neutral fraction xm = 
PHi/pH = Qgas/^b gives us *Hi = bOQ.gash'^{0.02/Qhh^) or 
xm ~ 2.45 X 10“^. We also assume Ts S> T 7 and hence we 
see the HI in emission. Using these we have 

PHi(k, z) = (1 -I- Ps-s{k, z) (19) 

where /3 = f /b. The fact that the neutral hydrogen is in 
discrete clouds makes a contribution which we do not include 
here. Another important effect not included here is that the 
fluctuations become non-linear at low 2 . Both these effects 
have been studied using simulations (Bharadwaj & Srikant 
2004). 

4 RESULTS 

We quantify the anisotropies in the HI clustering pattern by 
decomposing w{9^^,9±) into different spherical harmonics 

^ y w{9, p) Pi{jl)dy (20) 

where Pi (il) are the Legendre polynomials and we have used 
w{9, jl) = w( 6 i||, 9±) where 9 = -|- 9^ and fl = d\\/9. The 

fact that PHi(k, 2 ) depends only on even powers of /i ensures 
that all the odd harmonics will be zero. In addition to the 
monopole wo{9), we have calculated the quadrupole W 2 { 9 ) 
and the hexadecapole 1 x 4 ( 9 ) and we use the ratios 1 x 2 /wo 
and 1 x 4 /Wo to quantify the anisotropies. 

4.1 Reionization 

We have restricted our analysis to 2 = 10 which corresponds 
to 129 MHz and assume that half the hydrogen is ionized 
at this redshift ie. fv = 0.5 (Zaldarriaga , Furlanetto & 
Hernquist 2004). Recent investigations indicate that at this 
redshift the comoving size of the ionized bubbles will be of 
the order of a few Mpc (Furlanetto , Zaldarriaga & Hern¬ 
quist 2004a). Further, Furlanetto , Zaldarriaga & Hernquist 
(2004a) also show that the bias is expected to have a low 
value (near unity) for large bubble size R. We consider the 
LCDM model with R — 3 h~^Mpc and be = 1.5 as the fidu¬ 
cial model (Model A) for which Figures 3, 4 and 5 show 
Wo, W 2 /W 0 and W 4 /W 2 respectively. To study how the signal 
depends on various factors, we have varied the background 
cosmological model and the parameters R and be- Further, 
we have also considered the signal without the effect of the 
peculiar velocities and the Poisson term in order to explic¬ 
itly isolate the contribution from these terms. Table 1 shows 


Table 1. This shows the different combinations for which we 
have calculated the signal expected at the epoch of reionization. 
Here R is in /r“^Mpc. Further, P.V and P.F respectively indicate 
whether the peculiar velocities and the Poisson fluctuations have 
been included. 


Model 

Background cosmology 

R 

be 

P.V 

P.F 

A 

LCDM 

3 

1.5 

% 

% 

B 

VDE 

3 

1.5 

V 

V 

C 

LCDM 

5 

1.5 

V 

V 

D 

VDE 

5 

1.5 

V 

V 

E 

LCDM 

3 

1.5 

V 

X 

F 

LCDM 

3 

1.0 

V 

V 

G 

LCDM 

3 

0.0 

V 

V 

H 

LCDM 

3 

1.5 

X 

V 

I 

LCDM 

3 

0.0 

X 

V 



Figure 3. For the different models shown in Table 1, this shows 
the monopole of the signal expected at the epoch of reionization. 
Here 1 corresponds to 1.9h~^Mpc in the LCDM. 


the various combinations for which we have calculated the 
expected signal. 

One of the salient features of the monopole wo (Fig¬ 
ure 3) is that the signal is dominated by the Poisson fluc¬ 
tuations from the discrete ionized bubbles on small scales 
whereas it traces the dark matter on large scales. The an¬ 
gular scale where the transition from Poisson fluctuations 
to dark matter occurs depends on the background cosmo¬ 
logical model. Further, the Poisson contribution increases if 
the bubbles are larger (increasing R) because the number 
density of bubbles falls. The dependence on the background 
cosmological model can be attributed to the fact that the 
comoving distance corresponding to 2 = 10 differs by around 
10% in the LCDM and VDE models whereby the same ion¬ 
ized bubbles correspond to different angles in the two mod¬ 
els. It should be noted that though the overall amplitude of 
Wo also changes with the background cosmological model, 
it would not be possible to distinguish this from a change 
in the neutral fraction which would have the same effect. 
Increasing be increases the overall-amplitude of the signal at 
large scales leaving it unchanged at very small scales. The 
effect of peculiar velocities, we find, depends on the value 
of be, and it increases the signal when 6c = 0 whereas it 
reduces the signal when be = 1.5. An interesting situation 
occurs when 6c = 1 where the coefficient of Ps-s{k) in eq. 
(18) nearly cancels out and the signal is extremely small on 
large scales. 

Considering next the anisotropy (Figure 4 and 5) we 
find that the dominant feature is a bump at small scales 
caused by the Poisson fluctuations from the discrete ion- 
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Figure 4. For the different models shown in Table 1, this shows 
the quadrupole to monopole ratio of the signal expected at the 
epoch of reionization. Here 1 corresponds to 1.9/i“^Mpc in the 
LCDM. 



Figure 5. For the different models shown in Table 1, this shows 
the hexadecapole to monopole ratio of the signal expected at the 
epoch of reionization. Here 1 corresponds to 1.9/i“^Mpc in the 
LCDM. 


ized bubbles. The location of the bump is sensitive to the 
size of the bubbles. The nature of the anisotropy is signifi¬ 
cantly altered if the Poisson fluctuations are not taken into 
account. The bump in the anisotropy is a consequence of 
the AP effect, as can be inferred from the fact that it is 
not changed much if the peculiar velocities are not taken 
into account. Further, this feature is not very sensitive to 
the background cosmological model. The bias be of the ion¬ 
ized bubbles changes the amplitude of the bump. The large 
scale anisotropy is a combination of both the peculiar veloc¬ 
ities and the AP effect which make opposite contributions 
to the anisotropy. The large scale anisotropy is nearly con¬ 
stant at a value which depends on both be and R. Finally 
we note the fact that the signal is highly anisotropic with 
W 4 /W 0 > W 2 /W 0 > 1, and it is possible that some of the 
higher angular moments not considered here have large val¬ 
ues compared to wo- 


4.2 Post-reionization 

We restrict our analysis here to z = 3.37 which corresponds 
to frequency 325 MHz for the HI radiation. In addition to 
the possibility of different background cosmological models, 
there is only one free parameter b namely the bias of the HI 
relative to the dark matter. We consider the LCDM model 
with 6 = 1 as the fiducial model (Model J) for which Figures 
6, 7 and 8 shows wo, 102/100 and w^/wo respectively. The 
other models which we have considered are listed in Table 2. 

The monopole wo (Figure 6) traces the dark matter 



arc minute 

Figure 6. For the different models shown in Table 2, This shows 
the monopole for signal expected at 325 MHz. Here 1 corresponds 
to 1.4/i“'^Mpc in the LCDM. 





Figure 7. For the different models shown in Table 2, this shows 
the quadrupole to monopole ratio of the signal expected at 
325 MHz. Here 1 corresponds to 1.4h“^Mpc in the LCDM. 


distribution. Though the amplitude varies in the different 
models not much significance can be attached to this as such 
an effect can also arise from variations in the neutral frac¬ 
tion which is quite uncertain. The anisotropy in the clus¬ 
tering pattern is a combination of the peculiar velocities 
and the AP effect. This is deduced from the large change 
in the anisotropy when the peculiar velocity contribution is 
dropped. Both the effects are sensitive to the the background 
cosmological model at this redshift (Figures 1 and 2). The 
point to note is that changing the bias b has an effect which 
is very similar to that of changing the background cosmo¬ 
logical model, and it will be hard to use observations of the 
anisotropies to individually constrain any one of them. 


5 DISCUSSION AND CONCLUSIONS. 

The AP effect and the peculiar velocities both introduce 
anisotropies in the redshift space HI clustering pattern. 
These two sources of anisotropy are parameterized by two 
different functions p{z) (Figure 1) and f{z) (Figure 2) re- 
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Figure 8. For the different models shown in Table 2, this shows 
the hexadecapole to monopole ratio of the signal expected at 
325 MHz . Here 1 corresponds to 1.4/i“^Mpc in the LCDM. 


Table 2. This shows the different combinations for which we have 
calculated the signal expected at post-reionization. Here b is bias. 
Further, P.V indicates the peculiar velocities of HI. 


Model 

Background cosmology 

b 

P.V 

J 

LCDM 

1 

V 

K 

LCDM 

1.8 

V 

L 

VDE 

1 

V 

M 

LCDM 

1 

X 


spectively, both of which depend on the background cosmo¬ 
logical model. We have focused on two different background 
models. One of which has Q,mO ~ 0.3 and Hao = 0.7 and 
another VDE model has 17^0 = 0.4 and Hao = 0.6. The two 
models differ in that in one of the models the dark energy 
density remains a constant (LCDM) whereas in the other 
model (VDE) it increases with redshift and saturates at a 
value thrice the present value at 2 ~ 3. We find that f{z) 
varies with the cosmological model only at low redshifts and 
it saturates at f{z) = 1 at high redshifts where the universe 
is dark matter dominated in most cosmological models. Thus 
the anisotropies introduced by the peculiar velocities will be 
sensitive to the background cosmological model only at low 
redshifts, and it cannot tell us much about the background 
cosmological model at high redshifts. The function p(z) is 
sensitive to the background cosmological models at nearly 
all redshifts z > 0.1, and the anisotropies introduced by the 
AP effect hold the possibility of allowing us to probe the 
background cosmology at high redshifts. 

We have considered simple models for the HI distribu¬ 
tion at the epoch of reionization and in the post-reionization 
era. These models have a few parameters which quantify our 
ignorance about the HI distribution during these epochs. For 
the reionization era we have assumed that the total hydro¬ 
gen content traces the dark matter, and the reionization 
proceeds through spherical bubbles of ionized gas of comov¬ 
ing radius R. Further, the centers of these bubbles could be 
biased with respect to the underlying dark matter. The HI 
in the post reionization era is assumed to be in high column 
density clouds which could be biased with respect to the 
underlying dark matter. For both the epochs we hnd that 
the anisotropies in the HI clustering pattern are sensitive 


to both, the background cosmological model and the free 
parameters determining the HI distribution. Given the fact 
that very little is known a-priori about the high redshift HI 
distribution, we conclude that it will not be possible to use 
observations of the anisotropies alone to uniquely constrain 
the background cosmological model. It may be noted that 
our findings contradict earlier claims (Nusser 2004) where 
it was proposed that the anisotropies in the epoch of reion¬ 
ization HI signal could be used to constrain the background 
cosmological model. 

We Hnd that the anisotropy in the HI clustering is a 
combination of the contributions from both the AP effect 
and the peculiar velocities, and it cannot be attributed to 
anyone of them alone. Further, the anisotropy pattern at 
the epoch of reionization is very sensitive to the reionization 
model. In our simple model where the ionized bubbles are all 
of the same radius, we find a bump in the anisotropy pattern 
at the scales corresponding to the size of the bubbles. The 
bump will possibly be smeared in a more realistic situation 
where the bubbles will have a spread of sizes, but we may 
still expect a prominent feature at the scale corresponding 
to the characteristic bubble size. The anisotropy at large an¬ 
gular scales, we find, is sensitive to the bias of the centers of 
the ionized bubbles. Our analysis leads us to the conclusion 
that observations of the anisotropy in the HI clustering at 
the epoch of reionization will be a very powerful probe of the 
size, shape and distribution of the ionized regions the uncer¬ 
tainties in these being much more than the uncertainties in 
the background cosmological model. In the post-reionization 
era, the changes in the value of the bias parameter have 
a similar effect as changing the background cosmological 
model. These observations can be used to constrain the cos¬ 
mological model if the bias parameter can be determined by 
independent means like the bi-spectrum (Verde et al. 1998; 
Scoccimarro 2000). 

Finally we note that the temperature two-point cor¬ 
relation function considered here may not be the optimal 
statistical tool for detecting and quantifying the high red¬ 
shift HI signal. The angular power spectrum (Zaldarriaga 
, Furlanetto & Hernquist 2004) and the visibility correla¬ 
tions (Bharadwaj & Sethi 2001, Morales & Hewitt 2004 and 
Bharadwaj & Ali 2004b) have been proposed as optimal 
statistics for this purpose. In this paper we have chosen 
to study the temperature two-point correlation because of 
its similarity to the galaxy two-point correlation function 
where the anisotropy introduced by redshift space distor¬ 
tions is well understood and has received much attention in 
the literature. Though both the AP effect and the peculiar 
velocities are both included in Bharadwaj & Sethi (2001) 
and Bharadwaj & Ali (2004b), how to make the best use 
of these observations to constrain background cosmological 
models and models for the HI distribution still remains an 
open issue. 
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